function y=Theorem4_3(u,mu,sigma,a,beta,gamma)
y=((2*mu/sigma^2-1)*(a-gamma)^(2*mu/sigma^2-2)/u^(2*mu/sigma^2-1))*(((2*mu/sigma^2-1)*beta^(-2*mu/sigma^2)/(a-gamma)^(1-2*mu/sigma^2))/(1-(beta/(a-gamma))^(1-2*mu/sigma^2))+((1-2*mu/sigma^2)*(1-(beta/a)^(1-2*mu/sigma^2))*beta^(-2*mu/sigma^2)/((beta/(a-gamma))^(1-2*mu/sigma^2)))/(1-(beta/(a-gamma))^(1-2*mu/sigma^2))^2)+(1-(u/(a-gamma))^(1-2*mu/sigma^2))*(-((2*mu/sigma^2-1)*beta^(-2*mu/sigma^2)/(a-gamma)^(1-2*mu/sigma^2))*((2*mu/sigma^2-1)*(a-gamma)^(2*mu/sigma^2-2)/beta^(2*mu/sigma^2-1))/(1-(beta/(a-gamma))^(1-2*mu/sigma^2))^2+(((1-(beta/a)^(1-2*mu/sigma^2))*(1-2*mu/sigma^2)^2*beta^(-2*mu/sigma^2)*(a-gamma)^(2*mu/sigma^2-2))*(1-(beta/(a-gamma))^(1-2*mu/sigma^2))-2*((2*mu/sigma^2-1)*(a-gamma)^(2*mu/sigma^2-2)/beta^(2*mu/sigma^2-1))*((1-2*mu/sigma^2)*(1-(beta/a)^(1-2*mu/sigma^2))*beta^(-2*mu/sigma^2)/((beta/(a-gamma))^(1-2*mu/sigma^2))))/(1-(beta/(a-gamma))^(1-2*mu/sigma^2))^3);
% A=(1-(u/(a-gamma))^(1-2*mu/sigma^2));
% B=(1-(beta/(a-gamma))^(1-2*mu/sigma^2));
% C=((1-2*mu/sigma^2)*(1-(beta/a)^(1-2*mu/sigma^2))*beta^(-2*mu/sigma^2)/((beta/(a-gamma))^(1-2*mu/sigma^2)));
% A'=((2*mu/sigma^2-1)*(a-gamma)^(2*mu/sigma^2-2)/u^(2*mu/sigma^2-1));
% B'=((2*mu/sigma^2-1)*(a-gamma)^(2*mu/sigma^2-2)/beta^(2*mu/sigma^2-1));
% C'=((1-(beta/a)^(1-2*mu/sigma^2))*(1-2*mu/sigma^2)^2*beta^(-2*mu/sigma^2)*(a-gamma)^(2*mu/sigma^2-2));
% k=((2*mu/sigma^2-1)*beta^(-2*mu/sigma^2)/(a-gamma)^(1-2*mu/sigma^2));
%y=A'*(k/B+C/B^2)+A*(-k*B'/B^2+(C'*B-2*B'*C)/B^3);